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Abstract 

The procedures for studying a single narrow s-channel resonance or 
nearly degenerate resonances at a lepton collider, especially a muon 
collider, are discussed. In particular, we examine four methods for de- 
termining the parameters of a narrow s-channel resonance: scanning 
the resonance, measuring the convoluted cross section, measuring the 
Breit-Wigner area, and sitting on the resonance while varying the 
beam energy resolution. This latter procedure is new and appears to 
be potentially very powerful. Our focus is on computing the errors in 
resonance parameters resulting from uncertainty in the beam energy 
spread. Means for minimizing these errors are discussed. The discus- 
sion is applied to the examples of a light SM-Higgs, of the lightest 
pseudogoldstone boson of strong electroweak breaking, and of the two 
spin-1 resonances of the Degenerate BESS model (assuming that the 
beam energy spread is less than their mass splitting). We also ex- 
amine the most effective procedures for nearly degenerate resonances, 
and apply these to the case of Degenerate BESS resonances with mass 
splitting of order the beam energy spread. 



1 Introduction 



The possibility of analyzing narrow s-channel resonances is considered to 
be one of the most important strengths of muon colliders, at present under 
consideration by different laboratories Most of the content of this note is 
indeed addressed to physics of particular relevance to future muon colliders 
(for general reviews see [0). The importance of analyzing accurately an s- 
channel resonance at a lepton collider and the physical interest to be assigned 
to the extracted information are, of course, not new subjects. Much literature 
has in fact been devoted to the problem of the Z resonance shape at electron- 
positron colliders and to the implications of its accurate measurement. Many 
of the special issues that arise when studying a narrow s-channel resonance, 
such as a light SM Higgs boson, have also been considered [Q. 

However, these latter studies assumed that the beam energy profile is 
perfectly known. If a resonance is so narrow that its width is smaller than 
or comparable to the beam energy spread, uncertainty in the beam's energy 
profile can introduce substantial errors in the experimental determinations 
of the resonance parameters. One of the main goals of the present paper 
is to assess these errors. To this end, we shall present a detailed discussion 
of methods for studying a narrow s-channel resonance or systems of nearly 
degenerate resonances at a lepton collider, and consider applications to par- 
ticular cases. We carefully analyze the errors in the measured resonance 
parameters (such as the total width and partial widths) that arise from the 
uncertainty in the energy spread of the beam. Procedures for reducing this 
source of errors are studied and experimental observables with minimal sen- 
sitivity to beam energy spread are emphasized. Our study will assume that 
the beam energy spread is independently measured, as will normally be pos- 
sible. Expectations and procedures for this determination at a muon collider 
will be noted. 

Although the parameters of a narrow s-channel resonance can be deter- 
mined with minimal sensitivity to beam energy spread by measuring the 
total Breit-Wigner area, the peak total cross section and the cross sections 
in different final state channels, these measurements are not all easily per- 
formed with good statistical accuracy. Generally, a more effective method 
for determining resonance parameters is to scan the resonance using specific 
on- and off-resonance energy settings. In this paper, we also discuss the 
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possibly superior new method introduced in Ref. |Q in which one operates 
the colhder always with center of mass energy equal to the resonance mass 
but uses two different beam energy spreads: one smaller than the resonance 
width and one larger. All of these possibilities will be compared. The most 
immediate application is to the study of a light SM-Higgs, where, as already 
known p, ^, accurate measurements at a muon collider might make it 
possible to distinguish a standard model Higgs from the lightest Higgs of a 
supersymmetric model. We shall also analyze in detail resonant production 
of the lightest pseudogoldstone boson of dynamical electroweak symmetry 
breaking models. Following that, we consider two almost coincident vector 
resonances of the Degenerate BESS model when their mass splitting is 
larger than the beam energy spread. 

We shall then examine the case of two nearby resonances with mass split- 
ting much smaller than their average mass, assuming that the energy spread 
of the beam is of the same order as the mass splitting. Different possible 
procedures for reducing the errors in the determination of the physical pa- 
rameters will be examined. As an application we shall discuss such a situation 
in Degenerate BESS. 

In Section 2 we present the analysis for a narrow resonance and in Section 
3 the applications we have mentioned for this case. Section 4 gives the 
analysis for nearly degenerate resonances, and Section 5 a corresponding 
application. 



2 Analysis for a narrow resonance 

In this Section we will discuss several ways of determining the parameters of 
a narrow resonance at a lepton collider. 

For a given decay channel, an s-channel resonance can be described by 
three parameters: the mass, the total width and the peak cross section. 
These quantities can only be accurately measured if the beam parameters 
are very well known. Measurement of the mass requires precise knowledge 
of the absolute energy value. Measurement of the width and cross section 
requires excellent knowledge of the energy spread of the beam and the ability 
to determine the difference between two beam energy settings with great 
precision. 
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All of these beam parameters must be known with extraordinary accuracy 
in order to study a very narrow resonance. For example, consider a light SM 
Higgs boson. In [0] it was found that the beam energy -Ebeam must be known 
to better than 1 part in 10^ and the beam energy spread AE'beam must be 
smaller than 1 part in 10^ in order to scan the resonance and determine 
its width and other parameters. Below, we shall find that errors in these 
parameters due to uncertainty in Ai^beam will only be small compared to 
statistical errors associated with the measurements of a typical cross section 
(equivalent to a rate in a particular final state channel) if the error in the 
measurement of Aii^beam is smaller than ~ 1%. 

The independent measurement of Ai^beam can have both statistical and 
systematic error. Let us first consider the impact of statistical errors. We will 
argue that statistical uncertainties in AE'beam (and also -Ebeam itself) should 
not be significant. Consider measurement of some resonance parameter p 
through a series of observations over a large number (S) of spills. In a year 
of operation at a muon collider there will be something like S = 10^ spills. 
Each spill {i.e. each muon bunch) will contribute to the measurement of p, 
but the value of p interpreted from the cross section a observed for a given 
spill will be uncertain by an amount 6p. This uncertainty arises (a) because 
of the limited number of events Ni accumulated during the spill and (b) due 
to uncertainties in -Ebeam and Ai^beam for that spill. We have, using the short 
hand notation of Ef, for -Ebeam, 

Sp 

where c = fj^, a = and b = and AE^ are measured by 

looking at oscillations in the spectra of the muon decay products arising due 
to spin precession of the muons in the bunch during the course of a very 
large number of turns around the ring. These measurements are completely 
independent of the measured rate(s) in question. The fractional statistical 
errors ^ and for each spill are expected to be small as we shall sum- 
marize below. As a result, the contribution from the ^ term in Eq. (|1|) will 
normally be much larger than that from the 6Eb and SAEi, terms unless the 
coefficients a and b are very large. We will later learn that a and b will always 
be small under circumstances in which systematic errors in E^ and AE^ do 
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not badly distort the parameter measurement. To the extent that this is 
true, the statistical error 5p/p after a large number of spills, computed as 



5p 
P 



-1/2 

(2) 



will be dominated by the usual c/y/^~Ni term. Note that if this term were 
altogether absent, then ^ would be proportional to l/\/S times the per spill 
error from 6Eb and 6AEb. Thus, even if this latter were quite substantial, it 
would be strongly suppressed after accumulating data over 10^ spills. 

Since -Ebeam and Ai^beam are expected to vary somewhat from spill to spill, 
what is important is the statistical error with which they can be measured 
in a given spill. This has been discussed in [|], ^, |10|, |Tl| . Very roughly, the 
frequency of oscillation in the signal of secondary positrons from the muon 
decays (which signal is sensitive to the precession of the naturally-present 
polarization of the muons) determines -Ebeam and the decay with time of 
the oscillation signal amplitude determines Ai?beam/-Ebeam- At a resonance 
mass of 100 GeV, [] a beam energy spread of A_Ebeam/-Ebeam ~ 3 x 10~^ can 
be achieved by an appropriate machine design while maintaining adequate 
yearly luminosity (> 0.1 fb~^). Further, in this case, for each muon spill 
[i.e. each muon bunch), one can determine the beam energy itself to 1 part 
in 10'' (5 keV) and measure the actual magnitude of Ai?beam/-£'beam with an 
accuracy of roughly 1.67%. For Ai?beam/-£'beam ~ 10~^, as might be useful 
for a much broader resonance than the SM Higgs boson and as would allow 
for substantially larger yearly integrated luminosity (> 1 fb^^), -Ebeam can be 
measured to 2 parts in 10^ (100 keV) and Ai?beam/-Eboam can be measured 
with an accuracy of roughly 0.2%. As described above, this level of statistical 
error for -Ebcam and Ai?beam/-£'beam will result in only a tiny 6p/p error for a 
given resonance parameter unless the coefficients a and b of Eq. (0) are very 
large. Thus, in what follows, we will analyze the parameter errors introduced 
by an uncertainty in Ai^beamZ-^beam that is assumed to be systematic in 
nature. The level of such uncertainty is not well understood at the moment. 
The energy spread will be affected by possible time dependence of other 
beam parameters, such as emmittance, by time-dependent backgrounds to 
the precession measurements, and by other sources of depolarization. Absent 

^The precise figures given here are from 
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a detailed design for the machine and the polarimeter, |Tl| states that it is 
"quite certain" that the energy spread can be known with relative systematic 
error better than ~ 1%, even if Ai^beamZ-f^beam = 3 • 10^^. 

Given the precision with which the beam energy can be determined, it 
is clear that the mass of a resonance will be precisely known. We will not 
consider here the statistical errors in the measurements of the relevant cross 
sections, but these will later become important when choosing among the 
different procedures that we will discuss. Our main focus will be on the 
errors in the determination of branching ratios and total widths induced by 
a systematic uncertainty in the energy spread of the beam. Therefore, we will 
discuss the possibility of reducing such errors, and also discuss measurements 
which are independent of the energy spread. As already noted, these errors 
are expected to be important only when the width of the resonance is smaller 
than or comparable to the energy spread of the beam. 

We assume that the resonance is well described by a Breit-Wigner shape. 
For a resonance R of spin j produced in the s-channel and decaying into a 
given final state F, one has 

.-(E) = 4.(2, + 1) \e._m4Upt^ 

where M and F are the mass and the width of the resonance, and E = ^/s. 
We will work in the narrow width approximation and therefore neglect the 
running of the width. We will consider also the total production cross section 

Y(R — > f+f-)T 

We assume that in the absence of bremsstrahlung the lepton beams have 
a Gaussian beam energy spread specified by Ai?beam/-Ebeam = 0.01 /?(%), 
leading to a spread in the center of mass energy given by 

aE = E ~ 0.007 E . (5) 

In the absence of bremsstrahlung, the energy probability distribution f{E) 
can be written in the scaling form 

f(E) = ^J^^), (6) 
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where Eq is the central energy setting and 

Jf{E)dE = l. (7) 

Accurate measurements of a narrow resonance require an accurate knowledge 
of both aE and of the shape function g. We will explore the systematic 
uncertainties associated with errors in determining aE, assuming that the 
shape function g is known to be of a Gaussian form. 

After including bremsstrahlung (at an e^e^ collider, beamstrahlung must 
also be included) it is no longer possible to write the full energy distribution 
in a scaling form; bremsstrahlung has intrinsic knowledge of the mass of 
the lepton, which in leading order enters in the form \og{E /rrii). Since the 
bremsstrahlung distribution is completely known it may be convoluted with 
the scaling form of Eq. (|^) to obtain a final energy distribution which depends 
upon both aE and the (very accurately known) machine energy setting. Some 
useful figures illustrating the effects of bremsstrahlung at a muon collider are 
Figs. 30, 31 and 32 in Appendix A of [Q. One typically finds that the 
peak luminosity at the central beam energy, Eq, is reduced to about 60-80% 
of the value in the absence of bremsstrahlung, and that dC/dE falls below 
10~^ dC/ dE\^^-^^^ at ~ O.Q-Eq. We will discuss results at a muon collider 
both before and after including bremsstrahlung, and show that the errors 
introduced by uncertainty in the beam energy spread can be very adequately 
assessed without including bremsstrahlung. 

The measured cross sections are obtained from the convolution of the 
theoretical cross sections [see (|^)] with the energy distribution: 

a^{E) = 47r(2j + l)T{R ^ tr)T{R ^ F)h{T, aE, E) (8) 

where 

From these equations we can immediately draw two useful consequences: 
• the ratio of a^ to the production cross section is given by 

F 

^ = Bp, Bp = B{R F) (10) 
and does not depend on a^; 
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• due to the normalization condition the integral of is independent 
of erg 

J a^{E)dE = J a^{E)dE. (11) 

Of course, one must be very certain that the integral over E covers all 
of the resonance and all of the beam spread, including the low-energy 
bremsstrahlung tail. 

Let us now consider the narrow resonance approximation. In this case, we 
can write ^ 

MV^(i?) = 7r(2j + l)B,n-BF ,2 2fA (12) 

_|_ ^^/4 

where we have scaled the energy and the total width in terms of the energy 
spread evaluated at the peak (here we take aE ~ ctm) 

(13) 

The convolution can also be written in terms of these variables, with the 
result (see Eq. (|)) 

2 

MVf (E) = 7r(2j + l)B,n-BF I g{x - x) / 2// ^' 

J x' + 7^/4 

= B,+,-Bf^{x,^) (14) 

In particular, for the production cross section, we get 

MV(E) = 5,+,-<l>(x,7) (15) 

Some simple results that are valid in the narrow resonance approximation 
are the following: 

• If r ^ (Tm (and F is also large compared to bremsstrahlung energy 
spread), then 

A measurement of the peak cross section (summed over all final state 
modes) gives a direct measurement of B^+i-. In addition, the resonance 
shape can be scanned by using a series of different energy settings for 
the machine and F can then be determined. 
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• If r <^ (Tm, then 



(17) 



where we have assumed a Gaussian form for the beam energy spread. 
(We also neglected bremsstrahlung, but this could be easily included.) 
In this case, the peak cross section determines Bi+i- ■ T provided ctm 
is accurately known, ctjv/ could in turn be measured by changing E 
relative to M. However, statistical accuracy may be poor because cross 
sections are suppressed by T/gm in comparison to those obtained if 

r > (jm- 

When scanning a resonance, the above limits make it clear that the cross 
section is largest and that systematic errors associated with ctm are smallest 
when (Jm is as small compared to F as possible. However, for a sufficiently 
narrow resonance, the luminosity reduction associated with achieving ctm 
values significantly smaller than F is often so great that statistical uncertain- 
ties become large. Thus, we are often faced with a situation in which F is 
comparable to or only a few times larger than aM- (For example, this will 
be the case when studying a light SM Higgs boson at a muon collider.) It is 
in this situation that systematic errors in the resonance parameters deriving 
from uncertainty in ctjv/ can be enhanced. 

We now discuss four ways of extracting Bi+i-, B^+i-Bp in a specific final 
state channel and F from the data and assess the extent to which their deter- 
mination will be infiuenced by uncertainty associated with an independent 
measurement of ctm- 



Scan of the resonance - If F > ctm, one of the most direct ways to deter- 
mine the parameters of an s-channel resonance is through a three-point scan 
In this method one measures the production cross section dc at three 
different energies. However, the position of the peak is independent of the 
energy spread, and we do not expect that uncertainty in the latter will induce 
an error in the mass of the resonance (we have checked explicitly that this is 
indeed the case). Therefore we will assume in the following that the mass of 
the resonance is known with high accuracy, and, as a consequence, we will 
take into account only two scan points; one will be chosen at the peak and 
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Table 1: Errors on T and B induced by A(Tjv//o"a/ = +0.05 (upper lines) 
and A(Tm/o"m = —0.05 (lower lines) evaluated by choosing one observation 
at the peak and the other at E — M = kT, for k = 1, 2, 3. Bremsstrahlung 
is neglected. 





ir 


2T 


3r 


r/cTM 


AB/B 


Ar/r 


AB/B 


Ar/r 


AB/B 


Ar/r 


1 


+0.220 


-0.205 


+0.171 


-0.160 


+0.114 


-0.101 




-0.140 


+0.198 


-0.119 


+0.153 


-0.089 


+0.094 


2 


+0.065 


-0.075 


+0.043 


-0.037 


+0.035 


-0.023 




-0.053 


+0.072 


-0.039 


+0.035 


-0.033 


+0.023 


3 


+0.031 


-0.036 


+0.023 


-0.017 


+0.022 


-0.013 




-0.027 


+0.034 


-0.022 


+0.017 


-0.021 


+0.013 


4 


+0.018 


-0.020 


+0.016 


-0.011 


+0.015 


-0.009 




-0.017 


+0.019 


-0.015 


+0.010 


-0.014 


+0.009 



the other one off the peak. We will shortly discuss how the results depend 
on the position of this last point. The parameters {B, T) are then extracted 
by a two parameter fit. Here, B = Bg+i- if one is able to sum over all final 
state modes, and B = B^+^-Bp if one focuses on a particular final state. In 
principle this problem has a unique solution, by deconvoluting the observed 
cross section. However, the error in the determination of cta/ (the energy 
spread at the peak of the resonance) induces an error on the determination 
of the parameters of the resonance. Assuming that the measured value of 
(Tc at a given energy E is given by o"c(-E, -B, r, cta/) (here again we use the 
narrow resonance approximation to put ap ~ ctm), the changes in the values 
of B and r that result if ctm is shifted by an amount Actm can be evaluated 
through the equation 

a,(E, 5, r, au) = ac{E, B + AB,T + AT, au + Acxm) (18) 

or, from Eq. (P!^), 

$(a;,7) =$(x + Ax,7 + A7) fl + ^) (19) 
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where 



Ax = x{aM + A(Tm) — x^cm) ~ —x 

A7 = li^M + AfiM) - 7(aM) ^ ( — ^ ) (20) 

If we measure the cross section at two different energies, and assume that 
Acjm is the same at these two nearby energies (as would be typical of a 
systematic error), we can easily employ Eq. (|19D to determine the induced 
fractional shifts AF/F and AB/B. We fix one of the two points at the 
peak {x = 0), and we let the other one vary between E — M = T and 
E — M = 3F. The induced systematic errors are illustrated in a series of 
Tables. In Tables |l| and we assume a Gaussian energy distribution without 
including bremsstrahlung and tabulate errors induced by Ao"A//o"Af = 0.05 
and 0.01, respectively. As clear from Eq. (|19|), the fractional errors on F and 
B depend only on the E — M choice and on the reduced width 7, but not on 
B. 



Table 2: Same as Table |1] for AaM/o'M = +0.01 (upper lines) and AomIom = 
—0.01 (lower lines). 





IF 


2F 


3F 


r/cTA^ 


AB/B 


AF/F 


AB/B 


AF/F 


AB/B 


AF/F 


1 


+0.036 


-0.040 


+0.029 


-0.031 


+0.020 


-0.020 




-0.033 


+0.040 


-0.027 


+0.031 


-0.019 


+0.019 


2 


+0.012 


-0.015 


+0.008 


-0.007 


+0.007 


-0.004 




-0.011 


+0.015 


-0.008 


+0.007 


-0.007 


+0.004 


3 


+0.006 


-0.007 


+0.005 


-0.003 


+0.004 


-0.003 




-0.006 


+0.007 


-0.004 


+0.003 


-0.004 


+0.003 


4 


+0.004 


-0.004 


+0.003 


-0.002 


+0.003 


-0.002 




-0.003 


+0.004 


-0.003 


+0.002 


-0.003 


+0.002 



As expected from our earlier discussion, the fractional errors become 
smaller for larger T/aM- The fractional errors are also decreased by moving 
the second point of the scan away from the peak (see Tables and ||). On 
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the other hand the statistical errors for ac are increasing as 



$(0,7) 



(21) 



becomes smaller. (Due to presence of background, the decrease is not simply 



proportional to cJc^M) / ac{E).) By fixing the amount of luminosity that 
one can safely lose without incurring substantial statistical error in o"c(-E), 
and for a given ratio F/ctm, one can, by looking at Fig. ^ fix the second scan 
point. For instance, for F = 2(Tm, allowing for a loss in luminosity of 60%, 
we see that the point to be chosen is x ~ 2 corresponding to E k, M + T . By 
choosing the second point of the scan aX E = M + 2F, the resulting behaviour 
of the systematic errors vs. 7 is illustrated in Fig. ^. For a given uncertainty 
on the energy spread the errors decrease for increasing 7. Therefore, as noted 
earlier, one should employ the smallest value of ctm consistent with having 
luminosity large enough to give small statistical cxc errors. For instance, at 
7 = 1 the errors on the resonance parameters are of order 2.5 -^ 3.5 times 
Actm/ca/, whereas for 7 = 3 they are down to 0.3 4- 0.5 times 1\ouIom- 

Of course, a full study of the optimization of this procedure in a concrete 
case requires also taking into account the third scan point. In fact, one has 
the further freedom of moving the position of this point, with corresponding 
changes in the errors. One can show that by taking the third scan point 
in a symmetric position with respect to the other off-peak point, one gets 
the same errors shown in Tables |I| and |^ (in the absence of bremsstrahlung) . 
By taking an asymmetric configuration the results are in between the ones 
obtained with the corresponding symmetric configurations. For instance, by 
taking one point at the peak, one at ii^ = M + 2r and the third one at 
E = M - F, we find, for 7 = 1 and A(Tm/(Tm = 0.05, AF/F = -0.174 and 
A5/5 = +0.186. 

Inclusion of bremsstrahlung might complicate the picture developed above, 
since the induced errors could, in principle, depend upon the resonance mass 
and the value of i?, even when holding F/ctm and [E — M)/aM fixed. To 
determine if the effects of bremsstrahlung are important for determining sys- 
tematic errors coming from uncertainty in (Tm-, we begin by considering again 
^(^M/fTM = +0.05. We adopt M = 100 GeV and consider R = 0.003%, 
0.03%, and 0.1% at a muon collider. The resulting errors are given in Ta- 
ble ^ in comparison to the no-bremsstrahlung results of the upper lines of 
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Table 3: Errors for F and B induced by t^aul^M = +0.05 for no 
bremsstrahlung (first lines, from Table |I]), and after including bremsstrahlung 
for R = 0.003% (second lines), R = 0.03% (third lines) and R = 0.1% (fourth 
lines). The energy profile before bremsstrahlung is taken to be a Gaussian 
specified by the corresponding to a given R taking M = 100 GeV. 





ir 


2r 


3r 


r/ (tm 




A "P /"P 

Ai / i 


A D / D 
I\d I D 


A "P /P 

Ai / i 


A D / D 

Aid/ jd 


A P /P 

Ai / i 


1 


+0.22 


-0.21 


+0.17 


-0.16 


+0.11 


-0.10 




+0.22 


-0.21 


+0.18 


-0.16 


+0.12 


-0.11 




+0.22 


-0.21 


+0.18 


-0.16 


+0.12 


-0.11 




+0.22 


-0.21 


+0.18 


-0.16 


+0.12 


-0.11 


2 


+0.065 


-0.075 


+0.043 


-0.037 


+0.035 


-0.023 




+0.067 


-0.078 


+0.046 


-0.044 


+0.040 


-0.032 




+0.065 


-0.075 


+0.044 


-0.041 


+0.037 


-0.029 




+0.067 


-0.077 


+0.045 


-0.042 


+0.039 


-0.031 


3 


+0.031 


-0.036 


+0.023 


-0.017 


+0.022 


-0.013 




+0.031 


-0.035 


+0.024 


-0.020 


+0.023 


-0.017 




+0.031 


-0.037 


+0.024 


-0.020 


+0.023 


-0.017 




+0.031 


-0.038 


+0.024 


-0.021 


+0.022 


-0.016 


4 


+0.018 


-0.020 


+0.016 


-0.011 


+0.015 


-0.009 




+0.019 


-0.022 


+0.017 


-0.013 


+0.015 


-0.011 




+0.018 


-0.021 


+0.016 


-0.012 


+0.015 


-0.011 




+0.018 


-0.020 


+0.015 


-0.012 


+0.015 


-0.011 



Table |l|. The induced errors computed including bremsstrahlung are quite 
independent of R and (within the numerical errors of our programs) are es- 
sentially the same as those computed in the absence of bremsstrahlung. This 
can be further illustrated in the limit of very small Actm/cm- In this limit, 
one can simply compute the errors induced in B and F by using a linear 
expansion of Eq. ( plSf ) and solving the resulting matrix equation. The results 
for lS.BjB and AF/F expressed as a coefficient times I\omI(^m are given in 
Table ^ for four cases. The first case is that of no bremsstrahlung. (Note 
that for Actm/o'm = ±0.01 the errors predicted are quite close to those of 
Table |.) The other three cases are for M = 100 GeV and R = 0.003%, 



13 



Table 4: Errors for T and B expressed as AT /T = crAaM/o'&i and 
AB/B = cbAgm I o'm in the limit of small AaMl<^M- Same entry order- 
ing and bremsstrahlung assumptions as in Table ^. 





ir 


2r 


3r 


r/ om 


Cb 


Cr 


Cb 


Cr 


Cb 


Cr 


1 


3.4 


-4.0 


2.8 


-3.1 


2.0 


-2.0 




3.6 


-4.3 
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0.46 
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0.48 
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-0.36 


4 


0.35 


-0.40 


0.31 


-0.21 
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0.39 


-0.45 


0.30 


-0.24 


0.30 


-0.26 




0.39 


-0.43 


0.32 


-0.23 


0.29 


-0.23 




0.35 


-0.40 


0.29 


-0.24 


0.28 


-0.20 



0.03%, and 0.1%, including bremsstrahlung. Once again, we see that the 
error coefficients computed including bremsstrahlung do not depend much 
on R and are essentially the same (within the numerical errors of our pro- 
grams) as the error coefficients computed without including bremsstrahlung. 
Thus, in what follows we will discuss errors from Aaul^yM without includ- 
ing the effects of bremsstrahlung. Once a particular resonance is discovered, 
the effects of bremsstrahlung upon the results presented here can be easily 
incorporated to whatever precision is required. 

As regards e~*"e~ colliders, at which bremsstrahlung is a more substantial 
effect, a study (similar to that above for the \i~ collider) shows that the 
Gaussian approximation for studying systematic errors from Aou is again 
reasonable. Although beamstrahlung is also important at an e^e~ collider, 
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we did not study its effects. However, current e+e^ collider designs are such 
that the energy spreading from beamstrahlung is typically comparable to or 
smaller than that from bremsstrahlung. 

An interesting general question regarding the scan procedure is how to 
optimize the choice of (Tm relative to F so as to minimize the net statistical 
plus systematic error. Let us consider using k = 1 for the scan, assuming a 
resonance mass of M = 100 GeV and that Agm/o-m is small (< 0.02), so 
that the linear error expansion using coefficients cr and is valid. A very 
rough parameterization for the induced error in T when 0.2 < T/cm ^ 3 and 



The optimal choice of R for determining the parameters of a given resonance 
depends upon many factors: how the machine luminosity varies with R; 
the variation of AaM/o'M with R; the variation of the Ao"M-induced errors 
as a function of R; the magnitude of the resonance width (in particular, as 
compared to o"m); and the size of backgrounds in the important final states to 
which the resonance decays. At a muon collider, R in the range 0.1%-i-0.15% 
is the natural result and allows maximal luminosity. Increasing R above 
this range does not result in significant luminosity increase; in contrast, the 
luminosity declines rapidly as R is decreased below this range. A convenient 
parameterization for the luminosity of an = 100 GeV muon collider, valid 



The specific coefficient in Eq. (^3]) represents the most pessimistic estimate 
for the instantaneous luminosities. On occasion we shall also discuss results 
assuming a factor of 10 larger coefficient, which we shall refer to as the 
optimistic luminosity estimate. We do not know what to expect for the 
variation of the systematic error in o"m as a function of R. In order to 
understand how the optimization might work, we have adopted on a purely 

^The Lyear parameterization interpolates the three results of Table 5 in pH| taken from 



A(Tm/(Tm < 0.02 is 




(22) 



for 0.003% <R< 0.12%, is: Q 




(23) 



0. 
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adhoc basis a form that mimics the per-spill statistical error: Q 



( 



R 



) 



-0.57477 



2 ■ 10 



3 



(24) 



0.12% 



This form corresponds to a decreasing fractional systematic error in aM as R 
increases, as seems a hkely possibihty. Finally, it is useful to recall the value 
of (Tm as a function of R: 



For R = 0.003%, the smallest R that is likely to be achievable, these quanti- 
ties have the values Lyear = 0.1 fb""*^, AaM/o'M = 0.0167, and ctm ~ 2 MeV 
(for M = 100 GeV). The rough implications of these dependencies are as 
follows. 

• For a broad resonance, defined as one with F ^ O.OOIM, one should 
operate the muon collider at its natural R value of order 0.12%. The 
AcTM-induced errors will be very tiny, both because F/um is very large 
and because Actaj/cta/ should be small for such R. A precision scan of 
the resonance is readily possible in this case. Further, the resonance 
cross section is maximal for large T/aM and both background and signal 
rates vary slowly ClS db function of R. 

• For a resonance with F < O.OOIM, one will wish to reduce R until 
r ^ cta/. The primary reasons to avoid F < ctm are to avoid the 
large AcrA/-induced errors summarized above and to maintain adequate 
sensitivity to the resonance width. Keeping T/aM > 1, one can combine 
the above formulae with the predicted variation of statistical error as 
a function of R at fixed L to determine the optimal R choice. To 
illustrate, we will focus on the statistical error for F obtained using the 
three-point scan with measurements at ii^ = M, M ± 2crA/. If T/aM is 
large (typically > 7 -v- 10) the signal cross sections at the three scan 
points are roughly independent of R (but only for R values smaller 
than or not too much larger than that corresponding to o"a/). Since 
the background cross section is also independent of R, statistical errors 

•^This form for the statistical AaM/o'M error interpolates the results for the two R = 
0.003% and R = 0.1% cases given in Sec. 4.3.2 of using a power law form. 
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would not vary much with R if the total integrated luminosity is held 
fixed. More generally, one will find that, at fixed L, AF/F oc W with 
p > 0. (This is a very crude representation valid only for a limited 
range of R.) The exact rate of increase depends upon both F/cta/ and 
the background level. Such cases will be discussed shortly. However, 
even if p = 0, AF/F also varies as l/y/L, where the available L is given 
in Eq. (p^). Defining / = _R/0.003%, the net result is that we can write 
(AF/F)stat. = Cstat./^~°'^^^- Meanwhile, if we insert the expressions for 
crj\/ and Ao-m/ctm in the expression for the Acrjv/-induced AF/F, we 
find (AF/F)syst. = Csyst./^'''^^, where (for the k = 1 three-point scan) 
Csyst. ~ 0.06 (IJ^Y'^^^ if M = 100 GeV; we have normalized F to the 
approximate width of a 100 GeV SM Higgs boson. If p — 0.337 < 0, 
as would apply for p ~ 0, then, if the statistical and systematic errors 
are added in quadrature, the opposite signs of their / exponents means 
that there will be a minimum in the combined error as a function of 
/. If we use as a reference the value /o such that Cgtat. = Cgyst., then 
one finds that the optimal / is such that f / fo = R/ Rq ^ 0.7 when 
p ~ (requiring F/(Jm(-Ro) ^ !)• In other words, one should choose a 
value of R somewhat smaller than that value which would yield equal 
statistical and systematic errors. 

However, large T/aM cannot be achieved for the very narrow Higgs 
and pseudogoldstone bosons discussed later in the paper. Even for the 
larger and Mpo values of possible interest for this optimization 
discussion, we have, at best, T/aM{R = 0.003%) ~ 3 -i- 5. Starting 
from R = 0.003%, one finds that, at fixed total L, AF/F increases quite 
rapidly with increasing R. As two examples, V jaui^R = 0.003%) ~ 3 
(~ 5) for a SM-like Higgs boson with Mh ~ 140 GeV (~ 150 GeV). 
In these two cases, (AF/F)stat. increases by a factor ~ \/lO (~ -\/2) as 
R increases from 0.003% to 0.01%. Representing this increase using a 
power law (it is actually somewhat faster than a power law), we can 
roughly write (AF/F)stat. ~ cR^ at fixed L, with p ~ 1 (p ~ 0.3) for 
Mh = 140 GeV {Mh ~ 150 GeV). (Note that for R above the 0.003% 
to 0.01% range, the p values are much bigger.) From the analysis given 
in the preceding paragraph, we immediately see that if p > 0.337, the 
best statistical and systematic errors will both be achieved by taking R 
as small as possible. In the SM Higgs case, F/cta/ > 5, corresponding 
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to Mh > 150 GeV would be required before p < 0.337 and there 
could be some possible gain from increasing R. Similarly, in the case 
of the narrow pseudogoldstone boson P° it is only at the very largest 
Mpo = 200 GeV mass considered that p falls just slightly below 0.337. 
In practice, our a priori knowledge of F (from the initial scan needed to 
pin down the precise location of the H or P") will be too imprecise to 
allow for such optimization; one should plan on operating the machine 
at -R = 0.003% if initial information indicates that we are dealing with 
a Higgs or pseudogoldstone boson. 

Measuring the on-peak cross sections for different ctm - In this tech- 
nique one presumes that M is already quite well known and that one has 
in hand a rough idea of the size of F. (This is the likely result after the first 
rough scan used to locate the resonance.) One then operates the collider 
at E = M for two different values of aju (spending perhaps a year or two 
at each value). The results of Eqs. (|1^) and (0) show that if cr™"^ ^ F 

and crj^/^ > F, then cTc«/'^)/ac(a]{}^^) = . The systematic error 

in F is then given by (AF/F),yst. = Act^S-'^^/ctS''^. In practice, (rfr/^^M'' 
will be limited in size. If we define a]^^^"^^^ = Cj^^^Cm (the geometric 
mean value) and compute = crc(cr"™)/o-c(o-"/'') = $(0, 7min)/$(0, 7max) 
(where 7max,min = F/o"^'^'''™") as a function of F, our ability to measure F in 
this way for any given value of c^'^^/a™" is determined by the slope |s| of 
ln[F/(T^}'^*'''^^] plotted as a function of ln[rc]. The relevant plot from Ref. 
appears in Fig. For a known a^^^^^^"^, the |s| at any F/o"^}"*''''' gives the 
relation (AF/F)stat. = |s|(Arc/rc)stat., where /S.rc/rc is computed by combin- 
ing the fractional statistical errors for crc(o"™'^) and (Tcicrf^'^) in quadrature. 
The point at which the magnitude of the slope, |s|, is smallest indicates the 
point at which a given fractional statistical error in the cross section ratio 
will give the most accurate determination (as measured by fractional error) 
of F/cT^j'^*''^'. We observe that V/af^^^^^ ~ 2 3 gives the smallest |s| (and 
hence smallest statistical error), although |s| at F/a^f*"'^' ~ 1 is not that 
much larger. As expected, the larger cr^^/a^'^, the smaller |s| at any given 
V /a'^^^^^^. The systematic error in F due to uncertainties in a"^^ and cr]{}™ is 
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obtained from 

$(0,-4117) $(0, 



r+Ar 



-rinin _|_ A /T-min 



$ 0,-4r $ 



r+Ar 



(25) 



In most instances, the cxJJ}™ measurement is not very sensitive to the precise 
value of cr™/°, in which case the result is (Ar/r)syst. = Aaf/'^/a'^f''. More 
generally (Ar/r)syst. = AaM/o'M to the extent that Ao-m/ctm is the same at 
the different aM settings (as might apply for systematic errors of a certain 
type). From this we see one clear advantage of this technique: the systematic 
error in F due to systematic error in ctm does not increase with decreasing 
F/cr]JJ^^'™". If statistical errors for measuring the two (Tc's can be kept small, 
this is a clear advantage of the technique as compared to the scan technique 
in which small F/ctm leads to large AaM-induced systematic error even if 
statistics are excellent for all measurements. We will shortly discuss issues 
related to statistical errors. 

To continue our analysis of systematic errors, let us note that, in the limit 
of very high statistics for the cross section measurements {i.e. zero statistical 
error for Tc), precise values of both ^^^^ = F/crJ(}™ and 7max = ^/<^Tt^ 
obtained by the above procedure. Further, adc'^j^) and a^la'^f'^) are given 
by i?$(0,7min) and i?$(0,7max)5 respectively, where B stands for B^+i- or 
Bi+i-Bp, depending upon whether we are looking at the total rate or the 
rate in some particular final channel F. As a result, there is no systematic 
error in B from uncertainty in cta/, only statistical error associated with the 
number of events observed. This is another advantage of this technique. 

Of course small systematic errors are not important if statistical errors for 
the technique are not also small. We summarize the considerations [^]. For 
the SM Higgs and the P° we found that it is best to use o"™"^ corresponding 
to R = 0.003% and corresponding to R = 0.03%, so that af/'^'''^ corre- 
sponds to R = 0.003% X v^To. The parameterization for the variation of Lyear 
given in Eq. (|2|) implies that Lyear = 0.1 fb"^ (0.47 fb"^) for R = 0.003% 
(0.03%). If, for example, F/a^'^*''^' = 1, one finds (Xe«/°)/c7c(cT]S-^^) = 4.5, 
implying that the signal rate S'((jj\/) = ivyear(crA/)o"c(crA/) is nearly the same 
for a'^^^ as for a"^^^. However, the background rate B is proportional to L 
and thus B/S is a factor of 4.7 times larger at a^^"^ than at cr™'^. Conse- 
quently, the statistical error in the measurement of addM^^) '^iH be worse 
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than for (Jc{(^^m^) for the same S. For a given running time at a given 
CTjvf, one must compute the channel- by-channel S and B rates, compute the 
fractional error in (Jc{(yM) for each channel, and then combine all channels to 
get the net cTc(crA./) error. This must be done for aM 
One then computes the net and net cTc errors as: 



0"?}™ and (Tm = (y 



M 




+ 



1 -2 



A 
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max\ 
M ) 



crAcr 



M 



AcTeffT 
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-1/2 



(26) 



The ratio of running times at a™"^ vs. cr^t^ cannot be chosen so as to 
simultaneously minimize the net t^Ocjoc and /\r^/r^. The former is mini- 
mized by running only at o"™'^, while the latter is typically minimized for 
< 1. For the SM Higgs, a good compromise is to take 



' M 



M 



1. As demonstrated in the next section, it turns out that 
for both the P° and the SM Higgs boson, the ratio Lyear(i? = 0.03%) /Lyear(i? = 
0.003%) ~ 4.7 and the predicted cross sections and backgrounds are such 
that this technique is very competitive with the scan technique as regards 
statistical errors for F. 

Of course, for resonances (such as those of the Degenerate BESS model 
considered later) that have fairly large widths, the normal scan procedure 
can achieve superior results to the rc-ratio technique. This is because there 
will be little sacrifice in luminosity associated with choosing an R such that 
Vjou is substantially larger than 1. Measurements on the wings of the 
resonance will have nearly the same statistical accuracy as measurements at 
the resonance peak. 

Measurement of F using the ratio of the t^t~ final state and total 
cross sections - First of all let us discuss the measurement of absolute 
and relative branching ratios. This is possible by simply measuring the cross 
sections in different final states, holding the energy fixed: crf^/crf^ = Bp^jBp^ 
and cf/cTc = Bp. As apparent from Eq. (p!oD, such cross section ratios 
do not depend on gm- Therefore, we can measure branching ratios and 



''in the scan procedure, there is a similar difficulty. There, Bj S \s large for the off-peak 
measurements. 
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ratios of branching ratios with no systematic error induced by the energy 
spread. Of course, to measure ctc, we must be able to sum over all final 
states for which Bp is substantial, taking into account backgrounds and 
systematics associated with correctly determining the relative normalizations 
of different final states. This might not be easy, and could even be impossible 
if the resonance has some effectively invisible decays unless the branching 
ratio for invisible decays can be measured in a different experimental setting 
(in particular, one in which invisible decays can be effectively 'tagged' by 
producing the resonance in association with some other particle). 

Given Bi+£-, the value of cxc and Eq. (|1^) can then be used to determine 
r if (Tm is known. The error on F induced by uncertainty in ctm, given the 
absence of systematic error in Bi+i- , is obtained from Eq. (|T^) with AB = 0, 



If we measure Uc and af ^ at the resonance peak we have Ax = and 
Eq. (|27|) then requires A7 = 0; in turn. 



Therefore, the systematic error induced in T does not depend on 7. On the 
other hand, if we move away from the peak, the situation changes drastically, 
as illustrated in Fig. ^ which gives AF/F vs. 7 for Aaul<^M equal to 0.05. 
This figure can be trivially scaled for different values of AaM/o'M- We see 
that above 7 2 the error decreases as one moves further away from the 
resonance peak. For values below 7 «i 2 the situation is more complicated, 
but for each choice of the energy there is a zero in the error. | Therefore, by 
opportunely choosing the energy of the measurement one could try to work 
in a region where the systematic error is very much reduced. For instance, if 
F ^ 1.4(Jm, by taking E = M + T , one has AF = 0. Notice that the locations 
of the zeroes do not depend on Aau- 

There are two potentially severe difficulties associated with this technique. 
First, in many instances the decay mode has a small branching ratio. 
If the resulting event rate in the final state is not large, the statistical 
error for af^^ will be large. Statistics for E significantly different from M (as 

really precise computation of the locations of these zeroes would need to include 
the effects of bremsstrahlung. 




(27) 



A7 = 



AF _ A(Tm 



(28) 
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required to take advantage of the zeroes discussed above) would be the first 
to become problematical. Second, we have already noted that measurement 
of (Tc may be quite tricky. One must correctly normalize different final state 
channels relative to one another and assume that all final state channels 
with substantial branching ratio are visible. Regarding the latter, there is 
the alternative of measuring the branching ratio for invisible final states using 
other experimental situations/techniques and then making the appropriate 
correction to compute the full cJc given the contributions to J2f'^c ^^at can 
be measured in the s-channel setting. 

Of course, in some instances B^+i- will also be measured at other ma- 
chines or via other processes. Also in this case, one can immediately deter- 
mine r given a measurement of (Tc- Alternatively, if Bp is also known for 
some final state F (as well as Bg+g~), /{Bi+£-Bp) will give a determination 
of r with systematic uncertainty from Actm as described above. 

Measurement of the Breit-Wigner area - As apparent from Eq. ([Tl|), 
the energy integral of the cross section does not depend on ctm, and, in the 
narrow width approximation, is given by 

J ac{E) dE = 47r2(2j + 1) T • Be+i- (29) 

Therefore, even if the energy spread is very poorly known, so long as it 
does not change as one scans over the resonance one could still measure 
the integral of the cross section, which is proportional to the product F • 
Bi+i-, and obtain Bi+g- from the ratio of the peak cross sections crf^^ /cXc 
(or, possibly, in another experimental setting). Of course, as discussed with 
regard to the previous procedure, the statistical error for af.^^ might be 
large even when measured at the resonance peak. In addition, the event rate 
substantially off-resonance (as compared to both F and cta/) would be small 
(while typical backgrounds would be essentially constant). Thus, one can 
often only measure the cross section with good statistical accuracy over a 
limited portion of the full range. 

If F ^ o"j\/, and aM is known ahead of time (via spin-precession measure- 
ments or the like), this fact will quickly become apparent after measuring the 
cross section at a few energy settings with E — M ^ aM- A reliable value for 
F ■ Bi+i- can normally be obtained so long as statistics are good at the peak. 



The procedure is that based on Eq. (p^Bf). One will measure (Jc{E) for a set of 
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points perhaps out to — M ~ ±2r, and the remainder of the integral will 
be determined using the Breit-Wigner shape that would have been revealed 
by the measurements made. Sensitivity to au will be very minimal. 

If au ^ r, then the energy dependence of adE) is determined by ctm, 
see Eq. (|l^, or more generally by g of Eq. (|]) (bremsstrahlung should be 
included). Given a known value for (Tm and a set of measured points, the 
remainder of the integral can be computed. However, statistical errors are 
likely to be quite large because the cross section is suppressed by the ratio 

r/a-M- 

Clearly, the most difficult case is that in which the resonance is very 
narrow and the smallest achievable ctm value is such that F ~ ctm- Unless 
statistics remain good far off the peak, which is not likely, deconvolution of 
the effects of au and T is required. The simplest deconvolution procedures 
for known are the scan and ratio procedures outlined previously. Thus, 
if r ~ (Ta/, the three-point scan technique and the technique of varying 
while keeping E = M are the most efficient for determining the parameters of 
an s-channel resonance so long as systematic uncertainty in the energy spread 
of the beam is smaller than 1% -h 2%. (Of course, the appropriate strategy 
for exploring a resonance would be quite different if ctm is not independently 
measured with high statistical accuracy using the precession measurements.) 



3 Applications 



SM-Higgs boson - In this Section we will apply the first two methods of 
Section H to the study of a light SM-Higgs boson. Consider first the three- 
point scan. By using the results of Fig. |^ one can easily determine the 
errors induced by the uncertainty in the energy spread as a function of the 
Higgs mass. For the evaluation of the Higgs width as a function of the mass, 



we make use of the expressions given in Ref. |T2[. The decay channels we 
have considered are bb, t^t~ , cc, gg and WW*, ZZ* (one of the two final 
bosons being virtual). The resulting total width is given in Fig. |^. We have 
considered the interval 50 < MH^GeV) < 150, but recall that, actually, there 
is an experimental lower bound at 90% C.L. of about 90 GeV Below 
110 GeV, the width of the Higgs increases approximately linearly with the 
mass (aside from logarithmic effects due to the running of the quark masses) 
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which means that the ratio Th/ctm is approximately constant. By choosing 
R = 0.003% (see 0) we get Th/ctm ^ 1- From Fig. |, we see that the Acta/- 
induced errors in B and in Th are about 15% and 2.5% for Acm/o'm = 0.05 
and 0.01, respectively, for a = 2 scan (as appropriate, given that ~ ctm, 
for comparing to the E = M, E = M±2ctm scan summarized below that was 
used to estimate statistical errors). Above 110 GeV, the systematic errors 
decrease rapidly due to the fast increase of the width. Fig. ^ summarizes 
the AB/B and ATh/^h fractional systematic errors as a function of Mh for 
R = 0.003% and k = 2. Notice that, at least in the region up to 110 GeV, it 
will be mandatory to have values of R of the order 0.003%. For instance, if we 
take R = 0.01%, corresponding to Tn/aM ~ 0.3 (for 50 < Mj^(GeV) < 110), 
the fractional AdAf-induced errors in B and Th increase to about 14%-^ 16% 
for AaM/o'M = 0.01, as can be seen from the dashed line of Fig. 0. 

Table 5: Percentage errors (la) for acB{H —>■ bb,WW*, ZZ*) (extracted 
from channel rates) and Th for s-channel Higgs production at the muon col- 
lider assuming beam energy resolution of R = 0.003%. Results are presented 
for two integrated four-year luminosities: L = 4 fb~^ (L = 0.4 fb~^). An 
optimized three-point scan is employed using measurements at E = Mh, 
E = Mh + 2a M and E = Mh — 2a m, with luminosities of L/5, 2L/5 and 
2L/5, respectively. [For the cross section measurements, this is equivalent 
to L ~ 2 fb"^ (L = 0.2 fb"^) at the E = Mh peak]. This table is taken 
from Ref. |]5|. Efficiencies and cuts are those employed in [^. The effects of 
bremsstrahlung are included. 



Quantity 


Errors 


Mass (GeV) 


80 


Mz 


100 


110 


cJcBihh) 


0.8%(2.4%) 


7%(21%) 


1.3%(4%) 


1%(3%) 


acB{WW*) 






10% (32%) 


5%(15%) 


acB{ZZ*) 








62%(190%) 


Th 


3%(10%) 


25%(78%) 


10% (30%) 


5%(16%) 


Mass (GeV) 


120 


130 


140 


150 


acB{bb) 


1%(3%) 


1.5%(5%) 


3%(9%) 


9%(28%) 


acB{WW*) 


3% (10%) 


2.5%(8%) 


2.3%(7%) 


3% (9%) 


acB{ZZ*) 


16%(50%) 


10% (30%) 


8% (26%) 


11%(34%) 


Th 


5%(16%) 


6%(18%) 


9% (29%) 


34%(105%) 
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It is interesting to compare these systematic errors to the statistical errors. 
The analysis at a muon collider done in Ref. gives statistical errors for a 
three-point scan using scan points at E = M, E = M±2ctm and R = 0.003%, 
assuming L = 4 fb~^ or L = 0.4 fb~^ total accumulated luminosity (corre- 
sponding to 4 years of operation for optimistic or pessimistic, respectively, 
instantaneous luminosity). The results of that analysis are summarized in 
Table ^ Except for Mh ~ Mz, the L = 4 fb~^ statistical error for measuring 
the total width would be of order 3% ^ 10% when 80 < Mj^ < 140 GeV. 
Therefore, to avoid contaminating this high precision measurement error with 
systematic uncertainty from Actm we will certainly want to have AaM/crn ^ 
1%. If one adopts the L = 0.4 fb~^ luminosity assumption (the benchmark 
value of Refs. [|r^, [TT]]), and if ttih < 130 GeV and not near Mz, the statisti- 
cal measurement error for F/^ is in the 10% to 20% range. This means that 
A(Tm/c"m as large as 5% would be very undesirable. Finally, we re-emphasize 
the fact that performing the scan using larger R than R = 0.003% leads to 
larger statistical errors until Mh approaches the WW decay threshold and 
ThI(Tm{,R = 0.003%) > 5. For lower Mh, the R = 0.003% resuhs are the 
best that can be achieved despite the smaller luminosity at R = 0.003% as 
compared to higher R values. For example, the error in Th for a given lu- 
minosity using R = 0.01% can be read off from Fig. 13 of [Q. One finds 
that L{R = 0.01%) /L{R = 0.003%) = 20, 10, 2 is required in order that the 
Th statistical errors for R = 0.01% be equal to those for R = 0.003% at 
Mh = 130, 140, 150 GeV, respectively. Existing machine designs are such 
that Lyear(^ = 0.01%) / Ly^^,{R = 0.003%) = 0.22 fb'VO-l fb"^ = 2.2. Thus, 
increasing R would not improve the scan-procedure statistical errors until 
Mh > 150 GeV. In addition, Aaju-induced systematic errors always rise 
rapidly with increasing R. For Mh < 150 GeV, one should employ the 
smallest value of R possible. 

Let us now compare to the rc-ratio technique. We employ the same total 
of 4 years of operation as considered for the three-point scan, but always 
with E = Mh- As noted earlier, a good sharing of time is to devote two 
years to running at R = 0.003%, accumulating L = 2 fb"^ (L = 0.2 fb"^) 
for optimistic (pessimistic) instantaneous luminosity, and a second two years 
to running at i? = 0.03%, corresponding to [using the luminosity scaling law 
of Eq. (|2|)] L = 9.4 fb"^ (L = 0.94 fb"^) of accumulated luminosity. For 
R = 0.03%, Tn/crjf'^ ~ 0.1 and, as noted earlier, each acBp cross section is 
decreased by almost the same factor by which the luminosity has increased. 
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Table 6: Percentage errors (la) for acB{H bb,WW*, ZZ*) (extracted 
from channel rates) and Th for s-channel Higgs production at the muon col- 
lider. Operation at E = Mh is assumed and the rc-ratio technique is used for 
determining Th- Results are presented assuming accumulated luminosities of 
L = 2 fb"^ (L = 0.2 fb"^) at R = 0.003% and L = 9.4 fb"^ (L = 0.94 fb"^) 
a.t R = 0.03%, corresponding to roughly two years of running at each R for 
optimistic (pessimistic) instantaneous luminosity assumptions. Efficiencies 
and cuts employed appear in [||. The effects of bremsstrahlung are included. 



Quantity 


Errors 


Mass (GeV) 


80 


Mz 


100 


110 


(JcBibb) 


0.7%(2.2%) 


6.3%(20%) 


1.2%(3.8%) 


0.9%(2.8%) 


acB{WW*) 






8.2%(26%) 


3.8%(12%) 


cTcBiZZ*) 








60%(190%) 


Th 


6%(19%) 


63%(200%) 


14%(45%) 


8%(25%) 


Mass (GeV) 


120 


130 


140 


150 


acB{bb) 


0.9%(2.8%) 


1.4%(4.4%) 


2.4%(7.6%) 


6.6%(21%) 


acBiWW) 


2.4%(7.7%) 


1.8%(5.7%) 


1.6%(5.0%) 


1.8%(5.6%) 


(JcB{ZZ*) 


15%(46%) 


7.9%(25%) 


6.3%(20%) 


7.0%(22%) 


Th 


6.3%(20%) 


6%(19%) 


5.4%(17%) 


4. 7% (18%) 



leaving the number of signal events unchanged. However, the background 
is increased by a factor of 4.7. A complete calculation is required. This 
was performed in [Q. The various acBp statistical errors are summarized in 
Table |^ along with the ATh/T h statistical error computed by combining all 
the listed final state channels and following the procedure of Eq. (26). We 
observe that the ratio technique becomes superior to the scan technique for 
the larger Mh values {Mh > 130 GeV). This is correlated with the fact that 
Th/cm™ (where cr]JJ™ is that for R = 0.003%) becomes substantially larger 
than 1 for such Mh- In particular, for larger Mh, T h / crft^"^^^ is in a range 
such that \s\ and, consequently, the error in Th will be minimal. Thus, the 
two techniques are actually quite complementary — by employing the best of 
the two procedures, a very reasonable determination of Th and very precise 
determinations of the larger channel rates will be possible for all Mh below 
2Mw- 

Finally, we again note that the r^-ratio technique has the advantage that 
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the AcTM-induced systematic error in Th is equal to AaM/crM and will, there- 
fore, be smaller by a factor of about 2.5 (if Mh < 120 GeV) for the rc-ratio 
technique than for the scan technique and that there is no AaAf-induced 
systematic error in the Bi+i-Bp determinations. Thus, although the statis- 
tical errors for the r^-ratio procedure are larger than for the scan procedure 
when Mh < 130 GeV, if the fractional error in aM is substantially larger 
than 0.01, the r^-ratio could have net overall (statistical plus systematic) er- 
ror that is smaller than the scan technique down to Mh values significantly 
below 130 GeV. 



The lightest PNGB - The s-channel production of the lightest neutral 
pseudo-Nambu-Goldstone boson (PNGB) (P°), present in models of dynam- 
ical breaking of the electroweak symmetry which have a chiral symmetry 



larger than SU{2) x SU{2), has recently been explored |TJ, |T5|. In the 
broad class of models considered in [Q, the is of particular interest be- 
cause it contains only down-type techniquarks (and charged technileptons) 
and thus has a mass scale that is most naturally set by the mass of the 
6-quark. Other color-singlet PNGB's will have masses most naturally set 
by mt, while color non-singlet PNGB's will generally be even heavier. The 
Mpo mass range, that is typically suggested by technicolor models [0, is 
10 GeV < Mpo < 200 GeV. 

Discovery of the P° in the gg ^ 77 mode at the Tevatron Run II 

and at the LHC will almost certainly be possible unless its mass is either very 
small (< 30 GeV?) or very large (> 200 GeV?), where the question marks 
are related to uncertainties in backgrounds in the inclusive 77 channel. Run 
I data at Tevatron can already be used to exclude a P° in the 50 — 200 GeV 
mass range for a number of technicolors Npc > 12 — 16. In contrast, an e^e~ 
collider, while able to discover the P^ via e'^e^ —>■ '~fP^, so long as Mpo is 
not close to Mz, is unlikely (unless the TESLA 500 fb~^ per year option is 
built or Ntc is very large) to be able to determine the rates for individual 
7F final states {F = bb, t^t~ ,gg being the dominant P° decay modes) with 
sufficient accuracy such as to yield more than very rough indications about 
the parameters of the technicolor model. The 77 collider mode of operation 
at an e~^e~ collider would allow one to discover and study the P^ with greater 
precision. 

A fi~^fJ^~ collider would play a very special role with regard to determining 
key properties of the P°. In particular, the P^, being comprised of DD and 
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EE techniquarks, will naturally have couplings to the down-type quarks and 
charged leptons of the SM. Thus, s-channel production — >■ P°) is 
predicted to have a substantial rate for ^/s ~ Mpo. Because the P° has a 
very narrow width (see Fig. |^), not much larger than that of a SM-Higgs 
boson of the same mass, in order to maximize this rate it is important that 
one operates the collider so as to have extremely small beam energy 

spread, R = 0.003%. For such an R, the resolution in = ^/s of the muon 
collider, aE, is of order cr^; ~ 1 MeV(-E/50 GeV), whereas the P° width, Fpo, 
varies from 2 MeV to 20 MeV as Mpo ranges from 50 GeV up to 200 GeV 
(small differences with respect to come from running fermion masses). 
Thus, aE < Fpo is possible and leads to very high P° production rates for 
typical — > P° coupling strength. 

Assuming that the P° is discovered at the Tevatron, the LHC or (as might 
be the only possibility if Mpo is very small) at an e~^e~ collider (possibly op- 
erating in the 77 collider mode), the fi'^fi~ collider could quickly (in less than 
a year) scan the mass range indicated by the previous discovery (for the ex- 
pected uncertainty in the mass determination) and center on ^/s ~ Mpo to 
within < (Tm- A first very rough estimate of Fpo would also emerge from 
this initial scan. One would then proceed with a dedicated study of the 
P^. One technique would be to use the optimal three-point scan fl^ of the 

resonance (with measurements at £^ = Mpo and E = Mpo ± 2aM using 
R = 0.003%). The three-point scan would determine with high statistical 
precision all the fi'^fi'^ —^P^-^F channel rates and give a reasonably ac- 
curate measurement of the total width Fpo. For the particular technicolor 
model parameters analysed in |]l4l, 4 years of the pessimistic yearly lumi- 



nosity (Lyear = 0.1 fb ^) dcvotcd to the scan yields the results presented in 
Fig. 19 of [|14|. Sample statistical errors for acB{P^ all) and Fpo taken 
from this figure at Mpo = 60 GeV, 80 GeV, M^, 110 GeV, 150 GeV and 
200 GeV are given in Table ^. 

In the analysis performed in [|T^, we did not consider the errors induced 
by a systematic error in the energy spread. In Fig. M we show the AaM- 



^This figure gives tlie errors before taking into account the possible variation of the 
luminosity with Mpa. If the muon collider is built so that the y/s =100 GeV luminosity 
is maximized, then Lycar will be smaller (larger) than the ^/s = 100 GeV value for smaller 
(larger) y/s. The relevant luminosity scaling laws are those given in Eq. (7.2) of Ref. [l4| . 
The effects upon the statistical errors quoted below of such luminosity scaling are given 
in Fig. 20 of pj], but will not be included in our discussion here. 
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Table 7: Fractional statistical errors (Icr) for acB{P^ all) (combining bb, 
T~^T~ , cc and gg tagged-channel rates — see [0) and Tpo for s-channel P° 
production at the muon collider. We compare results for an i? = 0.003% 
three-point scan with total integrated luminosity of L = 0.4 fb^^ (corre- 
sponding to four years of running at pessimistic luminosity, with distribution 
L/haX E = Mpo, 2L/5 at E = Mpo + 2a m and 2L/5 at E = Mpo - 2a m) to 
results obtained using the rc-ratio technique assuming accumulated luminosi- 
ties at E = Mpo of L = 0.2/ fb"^ at R = 0.003% and L = 0.94(2 - /) fb"^ 
at -R = 0.03% (corresponding to roughly 2/ years of running at i? = 0.003% 
and (4 — 2/) years of running at R = 0.03% for pessimistic instantaneous 
luminosity assumptions). / (tabulated below) is chosen to minimize the error 
in Tpo. We employ the efficiencies, cuts and tagging procedures described in 
1^ . The effects of bremsstrahlung are included. 



Quantity 


Errors for the scan procedure 


Mass (GeV) 


60 


80 


Mz 


110 


150 


200 


CTcB 


0.0029 


0.0054 


0.043 


0.0093 


0.012 


0.018 


T po 


0.014 


0.029 


0.25 


0.042 


0.052 


0.10 


Quantity 


Errors for the rc-ratio procedure 


Mass (GeV) 


60 


80 


Mz 


110 


150 


200 


/ 


0.8 


0.7 


0.6 


0.8 


0.9 


1.0 




0.0029 


0.0062 


0.055 


0.010 


0.011 


0.016 


T po 


0.014 


0.028 


0.24 


0.041 


0.039 


0.053 



induced fractional errors for the branching ratio B (where B = Bi+^-Bp if 
we focus on a given final state or B = Bi+i- if we sum over all final states) 
and for Tpo predicted for the same three-point scan measurement {E = 
M, E = M±2aM, R = 0.003%) as employed for the statistical error analysis 
summarized above. Results are shown for the two cases of Aau/crM = 0.05 
(solid line) and AaM/o'M = 0.01 (dashed line). For R = 0.003%, we find 
AFpo/Fpo ~ cr (J Acta/ /o" A/ with cr □ ~ 1.5 for Mpo ~ 60 GeV falling to 
cr^o ~ 0.45 for Mpo ~ 200 GeV. For Mpo < 80 GeV, the result is that 
the induced AFpo/Fpo systematic errors are comparable to the expected 
statistical errors even for A(Tm/c"m = 0.01. For example, at Mpo = 60 GeV 
both the systematic error and the statistical error are of order 1.5%. In the 
neighborhood of the Z peak the errors from the optimal three-point scan 
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1^ are largely dominated by the Z background and the A(Tj\/ / ou effect can 
be neglected if /^au/cTM is not large. For Mpo ~ 150 GeV (~ 200 GeV), 
cr^o ~ 0.6 (0.45), i.e. AT po /T po ~ 0.6AaM/o-Af (0.45 Actm/ctm), while the 
statistical ATpo/Tpo error from the three-point scan would be of the 
order 5% (10%). Thus, for A(Tm/o"m < 0.01 the systematic error would be 
much smaller than the statistical error. As a result, for Mpo > 100 GeV, 
as far as the AaM-induced errors for Tpo are concerned, one could consider 
employing a value of R larger than 0.003%. As an example, R = 0.01% 
could be chosen. Since Tpo/aM is still > 1 for this R and Mpo > 100 GeV, 
we can expect that systematic errors will still be under control. The actual 
systematic errors resulting from an £^ = M, E = M ± 2(Jm three-point scan 
appear in Fig. g For Mpo ~ 150 GeV (~ 200 GeV) and AaM/^M = 0.01, 
one finds systematic error of AFpo/Fpo ~ 0.028 (0.022), which is indeed an 
acceptable level. However, as described in the previous Section, despite the 
factor of 2.2 increase [see Eq. (p3D] in yearly luminosity achieved by increasing 
R from 0.003% to 0.01%, the decrease in the signal to background ratio is very 
substantial and would lead to worse statistical errors unless Mpo > 200 GeV 
[for which Tpo/aM{R = 0.003%) > 5]. Thus, for Mpo < 200 GeV and typical 



model parameters as embodied in the choices of Ref. \T^\, one should employ 
R = 0.003% for the scan. 

Let us now consider the rc-ratio technique for the P°. We will compare 
to the scan technique using the choices R = 0.003% for cr]{}™ and R = 0.03% 
for cjJJ}''^. This means alf^'""^ ~ 6.3 MeV(Mpo/100 GeV). From Fig. ^ we 
then find Fpo/a^f*"'^' ~ 2/3 at Mpo = 50 GeV rising slowly to 1.6 at Mpo = 
200 GeV. This region is that for which the slope \s\ (see Fig. ^) is smallest. 
Consequently, the error in Fpo will be small if that for Tc is. The Actc/ o"c errors 
for 3 years of operation at pessimistic instantaneous luminosity (L = 0.3 fb~^) 



at -R = 0.003% were given in Fig. 19 of Ref. [|T^. We rescale these errors 
to L = 0.2/ fb^^ (corresponding to 2/ years of operation at i? = 0.003%), 
where / will be chosen to minimize the error in Tc- We also compute Aac/cc 
for L = 0.94(2 - /) fb"^ devoted to R = 0.03% running (corresponding to 
4 — 2/ years of operation at this latter R). The net Aadoc is computed using 
Eq. ( p6D after combining all final state channels. We also compute Atc/tc 
according to the Eq. ( p6|) procedure. The corresponding statistical error for 
Fpo is then computed using the appropriate \s\ slope value. We then search 
for the value of / (see above) such that AFpo/Fpo is smallest. The value of / 
and the corresponding errors for the combined-channel cXc and for Fpo error 
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are given in Table |^ for the same Mpo values as considered for the three- 
point scan procedure. For cTc, the rc-ratio procedure statistical errors are 
very similar to the 4- year three-point scan statistical errors for all Mpo values 
considered. For Fpo, the rc-ratio procedure statistical errors are as good as 
the tabulated 4-year three-point scan statistical errors for Mpo < 110 GeV, 
and become superior for larger Mpo values where Fpo/ a^™*'''^' is significantly 
larger than unity. 

As we have discussed, for the rc-ratio procedure, the fractional systematic 
error in Fpo is equal to that in ctm and there is no systematic error m. B = 
B^+^-Bp. In contrast, we have seen that for Mpo < 80 GeV the systematic 
errors in Fpo and B from the scan technique will be somewhat larger than 
AaM/o'M- Taking Fpo as an example, the scan systematic errors are of order 
I.SActm/ctm- As summarized earlier, this means that the scan systematic 
error for AcrAf/aA/ = 0.01 is essentially the same as the scan statistical error 
computed assuming pessimistic luminosity. For optimistic luminosity the 
scan systematic error would be dominant. Thus, the rc-ratio procedure will 
actually give better overall (systematic plus statistical) error for Fpo and, 
especially, B than the scan procedure for low as well as high Mpo values. 
This will become especially important if the luminosity available is better 
than the pessimistic value or if Actm/cta/ > 0.01. If a narrow resonance is 
observed in the 77 final state at the Tevatron or LHC, and if it has the weak 
coupling to ZZ that is typical of a pseudogoldstone boson, then the rc-ratio 
technique for precision measurements of its properties at the muon collider 
is strongly recommended. 

Degenerate BESS - The Degenerate BESS model describes two isotriplets 
of nearly degenerate vector resonances (L, R) characterized by two parame- 
ters {M,g"), the common mass (when the EW interactions are turned off) 
and their gauge coupling. The main feature of the model is the decoupling 
property, which implies very loose bounds from existing precision experi- 
ments [|13| as shown in Fig. ^. Another consequence is that the decay of 
the resonances into pairs of ordinary gauge vector bosons is quite depressed. 
The total widths for the two neutral states (L3, R3) are given by 

Fl3,p3 = M hLs,Rs {g/ g") , (30) 

where g is the weak coupling constant. The behaviour of the total widths as 
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functions of g/ g" is shown in Fig. |T0|, whereas for gjg" 1 one has 



/.i3 ^ 0.068 (^^^ , hn^^^^i^ . (31) 

The ratio of the widths, r^g/F/jg, in the interval < g/g" < .5 approxi- 
mately varies between 0.15 and 0.07. The weak interactions break the mass 
degeneracy giving rise to the mass splitting 

= Ml, -Mr, = M {g/g") . (32) 

The behaviour of for g / g" <^ 1 is 




The branching ratios into charged lepton pairs are almost parameter inde- 
pendent and rather sizeable 

5(1.3 ^ ^+^") ~ 4.5%, 5(i?3 ^ tt) ^ 13.6% . (34) 

In this Section, we will assume that um is much smaller than M*, and there- 
fore we can apply the previous analysis to the two resonances separately. 
The case au ~ M'^ ^ ^Lz,Rz will be studied in Section 4. By combining 



Fig. with Fig. |T0|, we easily obtain the AaM-induced fractional errors in 
the branching ratios and in the widths of the two resonances L3 and as 
functions of g/g" for different choices of ctm [equivalent ly, R] see Eq. (||)]. 



The results are given in Figs. 11 and 12 for L3 and i?3, respectively, for the 



choices A(TAf/o"A/ = 1% and 5%, and for R = 1% and 0.1%. (Because the res- 
onance widths are large, we do not need the very small values of R required 
to study the SM Higgs or the P°.) Combining these results with the bounds 
on the portion of parameter space still allowed by precision experiments, one 
can put lower limits on the masses of the resonances such that they have 
not been excluded and yet one is able to measure the widths with no more 
than a given systematic uncertainty from Actm- For instance, in the case of a 
machine with R = 1% (typical of an e~^e~ coUider) and for Actm/ctm = 1%, 
one finds, from Fig. [1^, that g/g" > 0.3 is needed to avoid Ao"A/-induced 
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errors in F^g above 2%. From Fig. the portion of parameter space that is 
still allowed by precision experiment can be roughly expressed by 



1000 GeV 



M 




\9' 



) 



(35) 



Therefore, the g/g" > 0.3 requirement converts to the requirement that 
M > 670 GeV. This means that, at lepton colliders with R = 1% and 
Ao-m/ctm = 1%, one will be able to measure F^g (for an allowed L3 res- 
onance) with a AcTAj-induced error of less than 2% only if the mass of the 
resonance is greater than 670 GeV. Similar considerations apply to R^, where 
keeping the systematic error in T^^ below 2% requires g/g" > 0.7, leading to 
M > 1540 GeV. In the case of a muon collider, R = 0.1% can be achieved 
while maintaining large luminosity. In this case, the AaM-induced fractional 
error for L3 is below 2%, for AaM/i^M = 1%, if g/g" ^ 0.12, which converts 
to M > 280 GeV for L3 resonances not already excluded by the precision 
data. The corresponding limits for R3 are g/g" > 0.3 and M > 670 GeV. In 
particular, we see from Fig. |^ that, at a machine working at the top threshold 
(about 350 GeV), only resonances with g/g" < 0.14 have not already been 
excluded by the precision data. Therefore, only a machine with R < 0.1% 
would be able to measure the width and branching ratios of an allowed reso- 
nance without encountering significant systematic uncertainty coming from 
Actm/ctm ~ 0.01. 



4 Analysis for nearly degenerate resonances 

In this Section we will discuss nearly degenerate resonances, i.e. the case 
of a mass splitting much smaller than the average mass of the resonances, 
= M2-M1 < M = {Mi + M2)/2. We will be interested in the case where 
the energy spread of the beam is of the same order of magnitude as the mass 
splitting between the two peaks, um ~ M*. We will also assume that the 
widths of the two resonances are much smaller than the mass splitting, i.e. 
Fi, F2 -C M*. It follows also that Fi, F2 <^ ctm- In this approximation, we 
can safely describe the cross section as the sum of two Breit-Wigner functions. 
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and furthermore we can use the narrow width approximation. Therefore, Q 

MVe = 5,V$(xi,7i) +5,\,-<l>(a;2,72) , (36) 

where Bl+^- and Bj+^- are the branching ratios of the resonances into , 
and 

xi = , X2 = = xi — m [67) 

and 

M2 = Ml + M^ m' = . (38) 

Cm 

We have also assumed M ^ Mi M2. For a Gaussian beam we recall that 
the function $(a;, 7) defined in Eq. (|^) is given by 

j(2j + l)/_ e-l-"/=^j^<i,. (39) 

We may evaluate this expression by performing the Fourier transforms of the 
Gaussian distribution and of the Breit-Wigner and then taking the inverse 
Fourier transform. We get 

$(x,7) = (2j + l)7r3/2^ dpe-P'^^+'P"" (e{p)e-^^'^ + 9{-p)e+^/^) . (40) 

Since we are assuming that the resonance widths are much smaller than 
the energy spread, we may evaluate Eq. (|40|) in the 7 — approximation, 
yielding 

$(x, 7) = (2j + l)7V27r=^/2g-xV2 ^ (4^) 

as earlier given in Eq. (p!7D for j = 0. As previously noted, this expres- 
sion shows that when the energy spread is much bigger than the width the 
convolution gives rise to a Gaussian function with spread ctm, and we loose 
any information about the total width. In fact, the total cross section Uc 
depends only on the product B^+i-'j = T{R —>■ i~^i~)/aM- In the case of two 
resonances with 71, 72 0, we thus find 



^We focus on the total cross section, but it should be kept in mind that we could also 
consider the cross section in a given final state. In this case, T{Ri £^£~) {i — 1,2) 
should be replaced by T{Ri £'^£~)Bp in all that follows. The AcrM-induced systematic 
errors on this product would be the same as for T(Ri — > £~^£~). 
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where 

T(Ri i+i-) 

= 7,5^+,„ = ^ . (43) 

Cm 



The function (^2]) is invariant under the substitution 

91 ^ 92, xi^m' - xi (44) 
The behaviour of the function is characterized by the ratio 

92 r(i?2 



= a (45) 



and by m'*. For small m**, the convolution of the Gaussian with the two 
Breit-Wigners has a single maximum in between and depending on the 
value of a. For instance, for a = 1 the maximum is at m'^ /2. In this situation, 
the second derivative of the function (^) has two zeros corresponding to the 
changes of curvature before and after the peak. By increasing the second 
derivative acquires a third zero (in fact, a double zero). This is due to the 
effect of the smaller Breit-Wigner which gives rise to a further change of the 
curvature. Just to get an idea, we list in Table 5, for several choices of a, the 
critical value ml of m'* at which this third zero occurs. As is increased 
further, there comes a point at which the two Breit-Wigner maxima start 
to show up. The minimum value of m'* required to see two maxima, mj, is 
given as a function of the ratio a in Table § Notice that the invariance (^) 
implies r?T.f(a) = m|(l/a). 

We can now discuss the type of measurements necessary to determine the 
parameters Mi, M2, gi and g2- In this discussion, we will use the notation 
X = Xi- We first determine the overall location in energy of the resonance 
structure by locating the absolute maximum of the cross section. For small 
m'*, such that the individual resonance peaks are unresolved, the cross section 
has a single maximum near the location of the resonance with the larger g; we 
will assume that it is g2 which is largest. For m'* large enough that the peaks 
are resolved, we may locate the larger maximum. We denote the location 
of the maximum in the cross section by -Emax, and choose this as our first 
energy setting. We write 

-E^max = CTAfa^max + Mi , (46) 



35 



Table 8: Values of at which the double zero of the second derivative of the 
function (|42|) occurs (mj(a)) and at which the two maxima of the function 
( ^21) start to show up (m|(a)). 



a 






1 


2 


2 


2 


1.85 


2.63 


3 


2.02 


2.85 


4 


2.20 


2.98 


5 


2.31 


3.08 



where Xmax is a function of a and which can be evaluated numerically 
from Eq. (^2]). For instance, for a > 2 it is a good approximation to assume 
s^max ~ (specially for bigger than the critical value m2), implying 
£^jnax — ^2- The value of the cross section at Xmax provides a second input 
into determining Mi 2 and gi^2- To complete the process of determining these 
four parameters, we need two more measurements. If is larger than m^, 
so that a second (lower) maximum is present, we can use the location of 
the second maximum and the cross section value at this second maximum 
as our two additional measurements. We effectively have four equations in 
four unknowns, two equations involving the derivatives of Eq. (|42|) and two 
involving the absolute magnitude of Eq. (^2]). If no second maximum is 
present, then we must effectively determine the slope of ac{x) of Eq. (|^) 
at some energy location away from Xmax and determine the cross section 
at this same location. Measurements of (Tc{x) at two nearby values of x 
away from Xmax are needed. That these approaches are really equivalent 
becomes apparent when we realize that the > m| procedure of locating 
the second maximum actually requires measuring the slope of o"c(x) and 
finding its second zero. (Recall that this discussion assumes infinite statistics 
so that we will end up effectively computing the minimum error that will be 
induced by uncertainty in ctm-) 

In the absence of a second maximum, the choice of the other two energies 
can be difficult if m'* is smaller than the critical value ml. In fact, in this case 
the convolution of the two Breit-Wigner looks very much like the convolution 
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with a single Breit-Wigner. Not surprisingly, to obtain reasonable errors it is 
necessary that cta/ be such that is at least bigger than m*. Let us assume 
that for m'* > ml we can approximately locate the energy corresponding to 
X ~ and that we measure the cross section at two points in its vicinity (as 
well as at X = Xmax)- If is greater than we may continue to employ 
the above procedure or we may use the alternative procedure outlined earlier 
based on the fact that the cross section has two maxima, one near x = and 
one near x = m^; in the alternative procedure, we measure the energy and 
the cross section at the two peaks. 

We consider first the former procedure that is the only choice if < 
We start our analysis by taking a = 2 and then later discuss the modifications 
for different values of a. From the measurement which fixes the energy scale 
(see Eq. (^61) ) we find the following relation between the parameter errors 
and the uncertainty in the energy spread: 

1 + ^^"j xmax(m^ + Am") + Ami - Xme^A^') = , (47) 

where AMi = aM^rrii is the error in Mi and Am" is the error in m". 
From this equation we can eliminate Ami in terms of the other errors. This 
has been done by using the approximation Xmax ~ The errors on M" 
and on the partial widths can then be determined using the three cross 
sections — crc(xmax) and <Jc{x) at two other x values near — following 
a procedure analogous to that discussed for a single resonance using Eq. 
(|18|). In particular, we assume measurements of the cross section at Xmax 
and at x = 0.1 and x = 0.2. The resulting fractional errors for T{Ri — > 
{i = 1,2) and M" are given in Fig. ^ as a function of M^/aM for 
Acm/o'm = 5%. As expected, the errors grow rapidly once m" falls below mi. 
As m" is increased above m", there is a change of curvature in the fractional 
error curves around the critical value m|, after which the fractional errors 
approach asymptotic limits. Notice that the asymptotic value of AM"/M" is 
zero, whereas AT{Ri i+i-)/T{Ri £+£-) AaM/crM for both i = 1,2. 
This is because the cross sections depend only on the ratio T{Ri aM- 



In Fig. |TJ we represent the same quantities but for a = 3, and we see results 
very similar to those for a = 2 except for changes due to the different values 
of the critical points m\ and mf^. We have tried different choices for the two 
measured points off the maximum, varying them up to x = 0.3, without any 
significant change in the results. 
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We now consider the alternative procedure outlined earlier that is pos- 
sible when two cross section maxima become visible, that is when the mass 
splitting is bigger than m|. In this case, the four measurements for deter- 
mining Ml 2 and gi^2 are the energy locations of the two maxima and the 
cross sections at these two maxima. The measurement of the energy of the 
second maximum gives rise to a condition similar to the one of Eq. p7|). The 
resulting fractional errors in ikf and T{Ri —>■ are given in Fig. |15] (for 

a = 2), and are essentially the same as obtained in the previous procedure 
when m** > m^. 

The basic conclusion from these analyses is that for of the order of the 
critical value m|, the fractional errors in the parameters of the resonances 
are of the order of the fractional error in ctm- For smaller the errors 
become very large. For significantly bigger than m|, the fractional error 
in the mass splitting rapidly approaches zero while the fractional errors for 
the T{Ri — i> partial widths become equal to the fractional error in 

the energy spread. In short, we can use the critical value m| in order to 
discriminate between a good and a bad determination of the mass splitting. 



5 Application to Degenerate BESS 

In Degenerate BESS one can show that the condition A/P ^ Tl, F/j is rather 
well satisfied (by one and two orders of magnitude respectively). Therefore 
we can apply the analysis of the previous Section. From Fig. |16|, we see that 
the value of a is almost constant and approximately 2.2 ioi g/g" up to 0.2, and 
then a increases up to ~ 4 for g/g" = 0.5. As discussed in the last Section, we 
can use the values of m2 given in Table |] in order to determine the minimum 
value of M^/aM needed in order to make a good determination of the mass 
splitting. For a = 2.2 one finds that the minimum value is = 2.68. As 
in our earlier single resonance discussions, for any fixed value of the energy 
resolution R, we convert this bound into lower bounds for g/g" and for the 
mass M of nearly degenerate resonances that have not already been excluded 
by precision experimental data. From M^/aM ^ 2.68 we get 

R=l%^^>0.16, M> 370 GeV , 
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R = 0.1% ^ 4 ^ 0.05, M ^ 130 GeV . 
9 

We see that for R = 1% a. machine with energy near the top threshold would 
just be at the border of being able to accurately measure the mass difference 
between the L3 and Rs resonances not already excluded by precision data. 



6 Conclusions 

We have considered the production of a narrow resonance via s-channel col- 
lisions of leptons {£ — e or n). Here, a 'narrow' resonance is defined as one 
that has width F substantially smaller than the beam energy spread AE'beam 
that is natural for the collider (and therefore is associated with the largest 
instantaneous luminosity). It will be convenient to use the parameterization 
A-E'beam = O.Oli? E'beam, where R is in per cent. For example, at a muon 
collider with center of mass energy E ~ 100 GeV, R ~ 0.12% allows for 
maximal L and L declines rapidly as R is forced to smaller values by com- 
pression techniques. A resonance with width F ^ O.OOIM would then be 
narrow. Our focus has been on the systematic error that might be intro- 
duced into measurements of the parameters of a narrow resonance due to 
systematic uncertainty in the value of AE^beamZ-Ebeam- The important pa- 
rameters that can be measured are the branching ratio of the resonance into 
the charged leptons {i.e. those that are being collided), the product of this 
leptonic branching ratio times that for the resonance to decay to a particular 
final state, and the total width of the resonance. We examined four methods 
for determining the resonance parameters: (1) a scan of the resonance; (2) 
sitting on the resonance and changing the beam energy resolution; (3) mea- 
surement of the cross section in the final state; and (4) measurement of 
the Breit-Wigner area. Methods (3) and (4) avoid the introduction of sys- 
tematic errors due to uncertainty in A£'beam/-E'beam, but for the integrated 
luminosities that are anticipated to be available the statistical errors asso- 
ciated with these techniques would be quite large for a narrow resonance. 
Methods (1) and (2) can provide resonance parameter determinations with 
small statistical error. However, even in the limit of infinite statistical accu- 
racy, determinations of the resonance parameters are sensitive to systematic 
uncertainties in AE'beam if T is not much larger than AE'beam- At a muon 
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collider, the smallest R that can be achieved is expected to be i? = 0.003%, 
for which Ai^beam ~ for a light SM Higgs boson and Ai^beam ~ T po/2 for 
the lightest pseudo goldstone boson of a technicolor model. Consequently, in 
these and other similar detailed assessment of the systematic errors 

in resonance parameter determinations introduced by uncertainty in Ai^beam 
is very important. 

We have performed a general analysis to determine the (systematic) er- 
rors in the measured resonance parameters induced by a systematic uncer- 
tainty in Ai^beam- We find that the induced fractional errors in the lep- 
tonic branching ratio (and also the product of the leptonic branching ratio 
times the branching ratio into any given final state) and in the total width 
can be expressed as universal functions of the ratio F/ctm, where ctm is the 
spread in total center of mass energy resulting from the beam energy spreads: 
aM/M = O.OIR/V2. In the case of the minimal three-point scan, with sam- 
pling at E = M, E = M ± kr, the error functions also depend on k. 

For a minimal three-point scan with k = 1, the induced AF/F fractional 
systematic errors were parameterized as a function of T/aM in Eq. (|2^). 
Very roughly, for T/aM ~ 2.5 we find that the induced fractional errors 
in F and B = Bi+£- or Bi+g-Bp (-F=final state) are of the order of the 
fractional uncertainty AaM/cTM- As T/aM increases above 2, the fractional 
errors smoothly decrease. For values of T/aM below 1, the fractional errors 
in the resonance parameters increase very rapidly. For example, for T/aM ^ 
1 (~ 0.2) the AcTAf-induced fractional systematic errors in the resonance 
parameters increase to ~ 3.5-f-4AcrA//o"A/ (~ 20-i-25AcrA.//o"A/), for the k = 1 
scan. Thus, to avoid large systematic errors from Actm, it is imperative to 
operate the collider with T/aM no smaller than 1. If R can be adjusted to 
achieve values significantly larger than 1, one can consider how to optimize 
the choice of R so as to minimize the total statistical plus systematic error. 
A discussion was presented leading to the following two basic conclusions, 
(a) For a broad resonance, defined as one with F ^ O.OOIM, one should 
operate the muon collider at its natural R value of order 0.12%. The AaM- 
induced errors will be very tiny, both because T/aM is very large and because 
AaM/cTM should be small for such R. (b) For a resonance with F < O.OOIM, 
one will typically wish to operate at an R that is significantly smaller than 
that value which would yield equal statistical and systematic errors. In a 
typical case, this would mean a value of R such that T/aM is larger than 
5 -7- 10. Of course, if the resonance is extremely narrow, it may happen that 
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T/ctm is of order, or not much larger than, unity even for R = 0.003%. In 
this case, it will normally be essential to run with R = 0.003% even though 
this R yields the smallest machine luminosity. Larger values of R lead to 
a drastic decline in the signal to background ratio in a typical final state 
that, in turn, leads to very poor statistical errors (given the rather slow 
compensating increase with R of the instantaneous luminosity). 

For a narrow resonance with F -C O.OOIM, the technique in which one 
sits on the resonance peak and measures the cross section for two differ- 
ent values of cta/ {cr^t"^ > F and cr™"^ < F) is a strong competitor to the 
scan technique. F/cx^'^*'''^^ (where a^;^"*'^^^ = \J cr'^^a^") is determined by the 
ratio Tc = ac{a'^™) / ac^cy^^"^), where Uc is the measured cross section. For 
'^m'^/'^m" order 5 to 20, the statistical error in F/cx^'^*'''^' is smallest if 
T/a^^^^^^ ~ 2 -h 3. The larger cr^^ / cr'^l'^ , the smaller the statistical error 
for F/cr^'^*'''^'. For a typical choice of CaT'^/^m™ = finds a statistical 

error of AF/F ~ LSAr^/rc for T /afl"^^^"^ ~ 2 ^ 3. This technique has the 
advantage that the AciM-induced systematic error in F is simply given by 
AF/F = A(Jm/c"m; while there is no systematic error in the determination 
of any of the B = B^+^-Bp branching ratio products (F = a particular final 
state). Of course, if the resonance is very narrow {e.g. as narrow as a light 
SM Higgs boson or a light pseudogoldstone boson), F/cr^'^*''^' ~ 2 3 will 
not be achievable. In this case, the best that one can do is to employ cr™" 
{(Ttr) as given by i? = 0.003% {R = 0.03%). The statistical error in F for 
such a situation is typically still very good. 

Let us now summarize how these results apply in the specific cases we 
explored. 

In the case of a three-point scan of the SM Higgs boson, we have shown 
that in the region Mh ^ 110 GeV it is mandatory to have R < 0.003%. In 
fact, even for this very small R value, Th/o^m is still < 1, the very minimum 
needed for accurate measurements of resonance parameters. For Th/c^m ~ 1 
the fractional systematic errors induced in Th from uncertainty in the beam 
energy spread are of order SAa^/aM for a A; = 2 scan. This should be 
compared to the typically-expected statistical errors tabulated in Table ^ 
For example, for Mh = 110 GeV the statistical error in ATh/^h is ~ 5% 
for optimistic 4-year integrated luminosity of L = 4 fb~^ at R = 0.003%; 
Actm/o-m ^ 0.01 would be needed for the systematic error to be smaller 
than the statistical error. For the pessimistic 4-year integrated luminosity of 
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L = 0.4 fb^^, the statistical error would be much larger (e.g. ATh/Th ~ 16% 
at Mh = 110 GeV) and the AdAf-induced error would be much smaller than 
the statistical error if Actm/o'm ^ 0.01. Note, however, that increasing R is 
not appropriate as this would push one into the Th/ctm < 1 region, implying 
large statistical errors and still larger systematic errors. 

For the cr^J}^'^'™™ on-peak ratio technique, one must choose cr]^™ corre- 
sponding to i? = 0.003% (TH/a^™ ~ 1). Results for statistical errors were 
presented in Table As a point of comparison, for optimistic (pessimistic) 
instantaneous luminosity and 4 years of operation, the net production rate 
error after summing over important channels is of order 0.8% (2.7%) for 
Mh = 110 GeV, and the ATh/Th statistical error is of order 8% (25%). 
Although the statistical AocjcTc fractional error is somewhat smaller for the 
ratio technique than for the scan technique, the IS^iijYii statistical error 
is larger. However, the ratio technique might still be better if AouI<^m 
were as large as 5%, especially if the optimistic luminosity level is available. 
This is because the systematic error in Vu is equal to Ao"a//(Tm for the ratio 
technique as opposed to "iAoulfyu for the scan technique. The rc-ratio tech- 
nique becomes increasingly superior as the assumed luminosity increases. For 
M-u > 130 GeV, the ratio technique gives smaller statistical errors for both 
(Tc-B and Yh than does the scan technique (for which statistical errors rapidly 
become very large). Indeed, the two procedures are nicely complementary 
in that at least one of them will allow a measurement of with statistical 
accuracy below 6% (20%) for optimistic (pessimistic) luminosity. 

For the lightest PNGB (P*^) of an extended technicolor model, the Agm- 
induced errors for the three-point scan method can be kept smaller than in the 
case of the SM Higgs boson. This is because, for typical model parameters, 
the P° has a width that is larger than that of a SM Higgs boson; Fpo/ gm > 2 



is quite likely for R ^ 0.003%. For example, the parameter choices of P 



give Fpo ~ 5 MeV vs. gm ~ 2 MeV at Mpo = 100 GeV. As we have seen, the 
resulting Ypo/gm ~ 2.5 yields AdM-induced resonance parameter fractional 
errors of order Agm/o'm- This can be compared to the statistical errors given 
in Table 0, computed assuming the pessimistic 4-year integrated luminosity 
of L = 0.4 fb~^. For example, at Mpo = 110 GeV the fractional statistical 
error for Fpo would be ~ 0.04, which is much larger than the systematic 
error if Agm/o'm ~ 0.01. For Mpo < 80 GeV, the statistical error in Fpo 
declines to the l%-i-3% level while the systematic error for Agm/^m = 0.01 
rises to about this same level. As Mpo increases from 150 GeV to 200 GeV, 
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the statistical error for Tpo rises from ~ 5% to ~ 10% while the systematic 
error is below 1% if AaM/o'M = 0.01. For optimistic integrated luminosity 
of L = 4 fb~^, the statistical errors would be smaller than quoted above. For 
AaM/o'M = 0.01, the induced systematic errors would generally dominate 
for Mpo < 80 GeV. 

The P° resonance is sufficiently narrow that one should also consider 
using the aJJJ^'''™'" on-peak ratio technique to determine Tpo. For 4- year 
pessimistic luminosity operation we find the statistical errors for Tpo given 
in Table |^. For the on-peak ratio technique, the systematic error in Fpo 
is equal to Aaul^M and for lower Mpo values would only be smaller than 
the statistical Fpo error if AgmIom 0.01. For optimistic luminosity, the 
1% systematic error induced in Fpo for AgmIom = 0.01 would dominate 
over the statistical error for all but Mpo ~ Mz and Mpo > 150 GeV. Most 
importantly, the statistical and systematic errors of the ratio technique are at 
least as good as, and often better than, obtained using the scan technique. 
For Mpo < 110 GeV, the statistical Fpo errors of the ratio technique are 
almost the same as obtained via the three-point scan (performed with R = 
0.003%), while the systemic Fpo errors from the ratio technique are smaller 
{Agm/o-m vs. ~ l.bAaM/oM)- For Mpo > 120 GeV, the Acrj\/-induced 
systematic errors are comparable for the two techniques, but the statistical 
errors for the ratio technique are substantially smaller than for the three- 
point scan. Precision measurements of the properties of a P° resonance 
would, thus, always be best performed using the ratio procedure. 

In the case of the two resonances of the Degenerate BESS model, we 
have determined (for several typical R values) the region of the model pa- 
rameter space for which the fractional errors in the resonance properties (F, 
. . .) induced by Agm are less than a given fixed value. Induced errors are 
small for large resonance masses. But, for any given choice of R, as the res- 
onance mass is decreased, while maintaining model parameter choices such 
that the model is still consistent with precision experimental data, there 
comes a point at which, even for AomI^m ~ 0.01, the induced error be- 
comes large. For allowed model parameter choices yielding a resonance mass 
below this, the resonance's properties cannot be measured accurately. The 
lowest masses of the resonances that correspond to precision-data-allowed 
BESS model parameters, and for which AcM-induced errors in the measured 
resonance properties are small, decrease rapidly with decreasing R. As a re- 
sult, the ability to achieve i? < 0.1% at a muon collider would be crucial for 
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exploring the low-resonance-mass portions of Degenerate BESS parameter 
space not currently excluded by precision data. 

Finally, we have performed the analysis of two nearly degenerate reso- 
nances, a situation encountered in a number of theoretical examples, includ- 
ing the Degenerate BESS model and the minimal supersymmetric model. 
We focused on the case in which the total widths of the resonances are 
much smaller than the mass splitting. We have shown that, in general, the 
AiTM-induced fractional error in the measured mass splitting, M*, and in the 
leptonic partial widths of the resonances (or leptonic partial widths times 
final state branching fraction) depend only on jou and on the ratio of 
the two partial widths, a. The main result is that the errors are generally 
big for I gm less than a certain critical value (typically in the 2 3 range) 
that is a function of a. As juu increases beyond the critical value, the 
AiTM-induced fractional error in the mass splitting approaches zero rapidly, 
whereas the fractional errors in the partial widths approach Ao"Af/<7M- As a 
concrete case, we have discussed the application to the two spin-1 resonances 
of the Degenerate BESS model for beam energy spread of the same order as 
the mass splitting between them. We determined the regions of the model pa- 
rameter space in which could be measured with AuM-induced fractional 
error below a given fixed value assuming a given value of R. The smaller the 
value of R that can be used while maintaining sufficient luminosity for small 
statistical errors, the larger the fraction of allowed parameter space for which 

can be measured with small AcrM-induced error. In particular, < 0.1% 
is required if we are to be able to separate the degenerate resonance peaks 
for model parameter choices not already excluded by precision experimental 
data in which the resonance masses are as low as 100 GeV. 
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Figure 1: Plot of the ratio ac{E) j ac{M) vs. x = {E — M)/aM, for T = kaM, 
with k = 0.5,1,2,3. 
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Figure 2: The fractional errors induced in B and F by Actm as functions of 
F/o"M when employing the three-point scan measurement {E = M, E = M± 
2F); results are shown for Acm/ctm = 0.05 (solid line) and AaM/c^M = 0.01 
(dashed line). 
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Figure 3: V/a'^'^'''^ plotted function of the cross section ra- 

ti o (Tc{(tZ"' )/(Tc{(Tm^) for various values of a'^^/aZ"' keeping a^''^''^^ = 
\I<^Z^<^ fixed. 
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Figure 4: The fractional error induced in F as a function of F/ctm for 
E — M = kr, with A; = 0, 1, 2, 3, assuming a known value of B^+i- from the 
peak measurements of cTc and af^^ . The value of AaM/o'M is fixed at 0.05. 
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Figure 5: The total width of the SM-Higgs boson as a function of Mh- 
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Figure 6: The AcxM-induced fractional errors in B and F as functions of Mh 
assuming a three-point resonance scan {E = M, E = M ± 2F); results are 
shown for AaM/o'M = 0.05 (solid line) and AaM/o'M = 0.01 (dashed line). 
The value of the beam resolution R has been fixed to 0.003%. 



53 



r(MeV) 



50 100 150 200 




50 100 150 200 



Mp^GeV) 



Figure 7: The total width of the hghtest PNGB P° as a function of Mpo 
(for the choice of model parameters given in Ref. [Q). 
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Figure 8: The AdAf-induced fractional errors in B and F = Fpo as functions 
of Mpo assuming a three-point resonance scan {E = M, E = M ± 2(7 m)'-, 
results are shown for f^ouloM = 0.05 (sohd line) and Acm/ctm = 0.01 
(dashed line), and for R = 0.003% and R = 0.01%. 
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Figure 9: The bound in the {M,g/g") parameter plane of the Degenerate 
BESS model arising from existing precision experimental data. The region 
to the left of the solid line is disallowed at > 90% CL. 



56 




Figure 10: The ratio F/M as a function of g/ g" for the Degenerate BESS 
for the neutral vector resonances L3 and R^. 
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Figure 11: The Degenerate BESS model: the fractional errors on B and F 
for the vector resonance L3 as functions of g/g" from the three-point scan 
measurement {E = M, E = M ± 2F), for A(Tm/o"a/ = 0.05 (solid line) and 
Ao-m/ctm = 0.01 (dashed line), and for R = 0.1% and R = 1%. 
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Figure 13: The fractional errors in the partial widths, Ar(i?j — > 
/V{Ri ^~^^~), i = 1,2 (dashed and dot-dashed lines respectively), 
and in the mass splitting, AM**/M* (continuous line), induced by an en- 
ergy spread uncertainty of Aum/cm = 0.05, as a function of for 
a = r(i?2 £~^i~)/T{Ri = 2. The input cross section mea- 



surements are at Xmax, see Eq. (Hq 



X = 
60 



0.1 and X = 0.2. 
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Figure 15: Same as in Fig. |T^, again with a = 2 but using the positions of the 
two maxima and the cross section values at the two maxima for determining 
the fractional errors. 
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Figure 16: The ratio r(L3 ^ £+r)/r(i?3 ^ vs. g/g" in the Degen- 

erate BESS model. 
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